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Abstract 

We show that exact results are obtained for the calculation of Casimir forces between arbitrary 
materials using the concept of surface impedances, obtaining in a trivial way the force in the limit 
of perfect conductors and also Lifshitz formula in the limit of semi-infinite media. As an example 
we present a full and rigorous calculation of the Casimir force between two metallic half-spaces 
described by a hydrodynamic nonlocal dielectric response. 
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In recent years experimental studies of Casimir vacuum forces between metallic surfaces 
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have reached the necessary accuracy to test in detail the theoretical predictions 
posed originally by Casimir in 1948 for perfectly-conducting parallel plates |5]. Indeed, the 
accuracy reached in the experiments has forced to consider detailed properties of the ma- 
terials such as absorptivity, rugosity, or finite temperature effects (f| Q, Isj]. The study of 
vacuum forces between realistic materials was pioneered by Lifshitz in 1956, who proposed 
a macroscopic theory for semi-infinite slabs described by a dielectric function e(uj) [10]. 
Lifshitz formula, which reduces to the Casimir result in the limiting case of perfect con- 
ductors, has been successfully employed in a number of experimental situations. Different 
authors have elaborated alternative derivations of Lifshitz formula that permit a simpler, 
more transparent approach to vacuum forces in realistic materials [llj|. Among several pro- 
posals to calculate the Casimir force the impedance approach was employed for the first time 
by Mostepanenko and Trunov |9| to derive the Lifshitz formula in an approximate fashion. 
In a series of investigations of vacuum forces in realistic materials, both at zero and finite 
te^erat™, KH-nchits^a Q et al. have found Sm all disomies ia the spatial he- 
haviour of the forces when calculated according to the impedance or the Lifshitz approach. 
In those papers it is argued that the discrepancies arise from the approximate nature of the 
impedance concept. Indeed, within the particular context of the optics of metals, the sur- 
face impedance is usually derived by taking approximations valid only for good conductors 
below their plasma frequency, and it is interpreted in terms of induced surface currents ^| . 
However, such approximations are unnecessary and the concept of surface impedance can 
be straightforwardly applied to arbitrary materials ^£,12, 3]. In an equivalent framework, 
Kats (l2j introduced the reflection amplitude coefficients r" (a = s,p) for metallic media to 
derive also an approximate version of Lifshitz formula to study the influence of nonlocality 
on van der Waals interactions in a semi-quantitative way. Noticeably, he stated incorrectly 
that the reflection coefficient cannot be expressed merely in terms of the surface impedance 
for dielectrics. Recently, a more rigorous derivation of the Casimir force in terms of the 
reflection amplitudes was discussed by Reynaud and collaborators |lj| using a S'-matrix 
formalism. 

The surface impedance Z is defined as the ratio of the complex electric and magnetic 



field components in the direction para 



wave, evaluated just inside the surface [!£ 



lei to the surface of incidence of an electromagnetic 



171 ] . For an s-polarized wave incident on a planar 
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surface z = zq moving (at an angle 9) towards the positive z direction, the impedance 
boundary condition is 

ECO = Z s (K t (z ) x z), (1) 
while for a p-polarized wave, the corresponding definition is 

Zm(zv) = z x E t (z ), (2) 

where z is a normal vector pointing inside the surface, and Zq denotes a position immedi- 
ately before the z = z interface. A main advantage of this concept is that it relates only 
external tangential fields, without the need of involving the internal degrees of freedom of 
electromagnetic field inside the material. Equations (0) and (J2J) are functional definitions of 
the surface impedance, so that they are exact and valid, not only for perfect conductors, but 
also for real metals and insulators. Indeed, from these definitions the reflection coefficients 
of the electromagnetic fields are obtained as 



Zt - Z s n _ Zl - zi 



o r p _ 



(3) 



Z° a + Z° a Z* + 2%' 

with Zq = q/ko, and Zq = k /q and Z? = k a c/(e a u), Z s a = u/k a c. These definitions in terms 
of the impedance are exact. For a local semi- infinite medium they yield the classical Fresnel 
coefficients 16^ . 

s ko k a p tako k a , , 

Ta ~ k^+K Ta ~ e a k + k a ' [ j 

However, the surface impedances are more general, as they are also valid for spatially dis- 
persive systems, for which the Fresnel relations are not applicable \v\ . Q] . 

In this paper we show that the use of Eqs. Q and (J2J) yields exact resuts for the Casimir 
forces using an expression for the force valid for arbitrary materials. As a non-trivial ap- 
plication of the surface impedance approach we calculate the Casimir force between two 
semi-infinite slabs with a non-local dielectric response. 

Consider two slabs a = 1,2 parallel to the xy plane within free space and separated 
by a vacuum cavity V of length L along the z-direction, with inner boundaries at z\ — 
and Z2 = L as shown in Fig.l. We assume that the slabs are non-chiral, translational 
invariant and isotropic within the xy plane, but otherwise they may be arbitrary. A given 
photon within V impinging upon a slab a may be reflected with a probability amplitude 
r" which depends on its polarization a, acquiring a phase kL as it moves on to the other 
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slab. Otherwise, it may be transmitted into the material with a probability T" = 1 — |r"| 2 
where it can be absorbed, exciting electronic or vibrational degrees of freedom, or it can be 
transmitted into the vacuum beneath the slab, in any case, becoming lost forever (multiple 
reflections within the slab are implictly accounted for in the reflection amplitudes r™. In 
thermodynamic equilibrium there would be photons coming from the outer vacuum and 
photons radiated by the materials themselves that woukd compensate exactly for the photons 
from v lost through absorption and transmission, appearing with a probability T* with no 
definite phase relation to the lost photons. Thus, the equilibrium radiation within the cavity 

V depends exclusively on its geometry, characterized by L, and on reflection amplitudes r". 
Thus, we construct an auxiliary system S' made of two infinitesimal sheets at z = z a , and we 
postulate that their reflection amplitudes are given exactly by the same amplitudes as those 
of the original slabs. This way, we assure that the radiation field within the fictitious cavity 

V corresponds to the real one. We further assume that in 5" a photon may be transmitted 
from V into the vacuum outside with an amplitude By choosing |r"| 2 + |t"| 2 = 1 we 
make certain that energy is conserved without having to account for any internal degrees of 
freedom within the fictitious sheets. 

We study first the case of s-polarized waves choosing the plane of incidence. 

Then the incident electric field can be written as E(r, t) = 'E e l ^ x ~ ult ^ 4>{z) , with amplitude 
E = (0, E y , 0), and the magnetic field B = (B x , 0, B z ), is determined by Maxwell equations: 
—iqB x = d z E y , and qB z = QE y , with wavevector q± = (Q, 0, ±k). The field component E y 
satisfies a one-dimensional wave equation. The solution satisfying the boundary conditions 
at slab 1 interface (z = + ) is 

E<(z) = e~ ikz + r{e ik \ (5) 
while at slab 2 interface [z = L~) we have 

E> (z) = e ik{z ~ L) + r s 2 e- ik{z ~ L) , (6) 

where the reflection amplitudes r" are given by Eqs.(OJ). 
The electric Green's function can now be calculated as 

G e M = SMSM, (7) 

where z< and z > are the smaller and larger of z and z', respectively, and W is their Wron- 
skian. The magnetic Green's function is obtained by replacing E y —>■ B x and r s a — ► — r s a in 



Eqs. ©-©• We do not consider B z separately, as it is simply prop ortional to E y . Therefore, 
for each Q, the local density of states per unit k 2 is given by yi| 

pU(z) = ~lm(G^(z,z)+Gf t (z,z)) , (8) 
with (k = k + i0 + ). By substituting Eqs. ©-© and its magnetic analogues we obtain 



p s k 2 = — ^Re 
2nk 



1 — r\r s 2 t 



(9) 



which is independent of z. Given the symmetries of the problem, the p-polarization density 
of states ff k 2 may be derived similarly, by just replacing B x — > —E x , and E y — > B y . The 
final expression is simply given by Eq.® after replacing all the superscripts s — > p. Finally, 
the total density of states is = p s k2 + p v k i- 

Consider now a photon in a state characterized by a, Q, and k 2 . Its momentum is ±hk 
and it moves with velocity ±ck/q along the z direction, so that its contribution to the 
momentum flux is hck 2 /q. Denoting by f(k) = N k + 1/2 to photon occupation number of 
state k, the momentum flux from the vacuum gap into slab 2 is then given by 

/ QdQdk 2 hc—f{k)p a k2 . (10) 

There is a similar contribution coming from the semi-infinite vacuum on the other side of 
the slab, obtained by substituting r% — > above and reversing the flux direction z —z. 
The total force per unit area is obtained by subtracting the contributions from either side: 



hc r dQQ f dk l Re i 



e - 1) + e 



(ii) 



A 2tt 2 J J q > q k 

with £° = (rfr^ exp(2ikL))^ 1 . The integral over k runs from iQ to and then to oo, so 
that q remains real and positive. This expression depends only on the reflection coefficients 
or, equivalently, on the surface impedances of the system and the slab separation. Lifshitz 
formula is recovered upon substitution of the local Fresnel amplitudes Eqs.(JU), whereas for 
perfect mirrors (r" = ±1) Eq. (jllj) yields the expected Casimir force. 

An example of a non-trivial application of the impedance approach is the calculation of 
Casimir force between two conductors with a dielectric function that shows spatial disper- 
sion. This is, e(r — r'jt — t'), or in the Fourier space e(q, u). Kats |3] studied nonlocal 



effects in an approximate way, and as he stated, it is necessary to specify correctly the de- 
pendence of the dielectric function on the wave vector. This is done in this section. In a 
conductor, the normal component of an incident p — polarized wave pushes the conduction 
charge towards or away of the surface creating an excess charge. When the frequency of 
the electromagnetic wave is greater than the plasma frequency u p of the metal, this charge 
accumulation propagates as a longitudinal wave (plasmon). Thus, at a surface, p— polarized 
waves couple to bulk plasmons taking away energy from the inicident wave and changing 
the values of the reflection amplitudes at the surface. In this case the medium is said to be 
active and it supports both longitudinal and transverse oscillations. We have to distinguish 
between the longitudinal dielectric function e(q, uS)i describing the response to the longitu- 
dinal dielectric field E|| (parallel to the wave vector) and the transverse dielectric function 
e(q, uj) t that is the response to a transverse electric field. The total field is the sum of the 
longitudinal and transverse parts E = E|| + Ej_. Thus, in the metal we have a p — polarized 
wave with wave vector (Q, 0, k) satisfying k 2 = e(u)tuj 2 /c 2 — k 2 and a longitudinal wave with 
wavevector 1 = (Q, 0, 1) that obeys the relation e(l,u)i = 0. Since s — polarized waves do 
not couple to plasmons the transverse dielectric function is local (no spatial dispersion). To 
describe the dielectric response of the metal we have a Drude transverse dielectric function 
et = 1 — uj 2 /{uj{ijj + ij)), where u> p is the plasma frequency and 7 the damping. For the 
longitudinal part we consider a hydrodynamic dielectric function [24] given by 

el ^ UJ ) = 1 -^—- ffl^.m . ( 12 ) 

ur + iwy — fj' L {Q z + l l ) 
where (3 2 = 3v 2 /5 with Vf the Fermi velocity of the metal 24j and it enters the hydrodynamic 
model as a compressibility term of the charge carriers an it is responsible for the spatial 
dispersion j^J. 

The surface impedance for p— polarized waves can now be calculated at the surface of the 
metal. However, since we now have the additional longitudinal field, we need an additional 
boundary condition. For the problem we are considering this condition is that the normal 
component of the total current be zero at the surface. This is, j z ± + j z \\ — 0. The surface 
impedance is calculated as 

E x ± + E x \\ 



Zp ~ B 



y 



kc | (e t -l)Q 2 c 
t t uo e t lio 



This expression for the impedance is now used in Eq.(j3J) to obtain the reflectivity for p — 
polarized waves in the non-local case as 

= £tfc — k — (e t — l)Q 2 /l , s 

Tp ~ e t k + k + (e t -l)QyV [ ' 

while the reflectivity for the s — polarized waves is the local Fresnel coefficient r s given in 
Eq.Q. In the local limit (Im(l) — > 0) the local Fresnel coefficient is recovered. 

The Casimir force between two semi-infinite conductors described by the HD model an 
now be calculated again using Eq. ((TTJ) . In Figure 3 we present the percent difference 
between the Casimir force for the local case(-Fk) with that obtained with nonlocal effects 
(F NL ), i.e. A% = \{F NL - F L )/F L \100. This is done for three different metals (K, Au, Al) 
as a function of the separation of the slabs in units of the plasma wavelength of the metal 
L/X p . The parameters for each metal (u> p , 7, Vf) are taken from the literature j^ll^j]. As 
expected the nonlocal effects become important for separations of the order of the plasma 
wavelength X p , and the difference between them can be significant for separations less than 
X p . For these separations the vacuum modes that contribute mainly to the Casimir force 
have a frequency larger than u p thus having propagating modes in the metals. Furthermore 
the better the conductor the higher the difference between the local and nonlocal cases. 

The non-coincident results for the Casimir forces associated to the Lifshitz and the 
impedance approach Refs. [l^ can be traced back to the introduction of an approximate 
expression for the impedance. In the original work of Mostepanenko and Trunov [9], they 



employed an expression for the impedance Z{uf) = l/ye^), valid only for small values of 
the electromagnetic wavevector Q, parallel to the surface of incidence. However, by intro- 
ducing the definitions Q and (J2J) in the same formalism, then the exact expression for the 
Casimir force fEq. (jlip ) is obtained by just replacing Z — ► Z p in Eq.(41) of Ref . ( ) . and 
Z — > Z s in Eq.(42) of the same reference. 

In conclusion, we have shown that with the correct definitions of surface impedance an 
exact result for the Casimir force is obtained. As case studies we calculated the surface 
impedance and Casimir force between an homogeneous slab with a non-local optical re- 
sponse. In this case, we show that nonlocal effects become relevant for separations less than 
the plasma wavelength of the slab. 
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FIG. 1: Schematics of the system used in our calculations. For the nonlocal case, the slab has an 
additional longitudinal field as shown in the Figure. 

FIG. 2: Percent difference between the Casimir force when local and nonlocal effects are considered. 
The separation between the slabs is in units of the plasma wavelength X p of the metals. The curves 
correspond to Au, Al and K. The horizontal line is a visual aid to show where the percent difference 
is 1%. 
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